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Abstract: In this work, the influence of gravity, fuel dilution, and inlet velocity on the structure, sta- 
bilization, and sooting behavior of laminar coflow methane-air diffusion flames was investigated both 
computationally and experimentally. A series of flames measured in the Structure and Liftoff in Com- 
bustion Experiment (SLICE) was assessed numerically under microgravity and normal gravity condi- 
tions with the fuel stream CH 4 mole fraction ranging from 0.4 to 1.0. Computationally, the MC-Smooth 
vorticity- velocity formulation of the governing equations was employed to describe the reactive gaseous 
mixture; the soot evolution process was considered as a classical aerosol dynamics problem and was 
represented by the sectional aerosol equations. Since each flame is axisymmetric, a two-dimensional 
computational domain was employed, where the grid on the axisymmetric domain was a nonuniform 
tensor product mesh. The governing equations and boundary conditions were discretized on the mesh 
by a nine-point finite difference stencil, with the convective terms approximated by a monotonic up- 
wind scheme and all other derivatives approximated by centered differences. The resulting set of fully 
coupled, strongly nonlinear equations was solved simultaneously using a damped, modified Newton’s 
method and a nested Bi-CGSTAB linear algebra solver. Experimentally, the flame shape, size, lift-off 
height, and soot temperature were determined by flame emission images recorded by a digital camera, 
and the soot volume fraction was quantified through an absolute light calibration using a thermocouple. 
For a broad spectrum of flames in microgravity and normal gravity, the computed and measured flame 
quantities (e.g., temperature profile, flame shape, lift-off height, and soot volume fraction) were first 
compared to assess the accuracy of the numerical model. After its validity was established, the influ- 
ence of gravity, fuel dilution, and inlet velocity on the structure, stabilization, and sooting tendency of 
laminar coflow methane-air diffusion flames was explored further by examining quantities derived from 
the computational results. 
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1 Introduction 

In 201 1, over 70% of electricity in U.S. was generated from the burning of fossil fuels, so a better 
understanding of combustion is of great importance [1]. Microgravity provides an ideal yet rigor- 
ous environment to study flame structure and examine computational models, because the structure 
of a microgravity flame is dramatically different from its normal gravity counterpart, and micro- 
gravity flames have larger scales and longer residence times, which makes it possible to study a 
wider range of flame conditions than under normal gravity. As a result, a significant amount of 
research has been conducted under microgravity conditions (e.g., [2-7]) to study the effects of 
various physical parameters on flames. In this study, previous computational and experimental 
investigations of coflow laminar diffusion flames (e.g., [2, 3, 8, 9]) were extended to explore fur- 
ther the influences of gravity, fuel dilution, and inlet velocity on the structure and stabilization of 
laminar coflow methane/air diffusion flames. 


2 Burner configuration 

The burner consists of a central jet, from which the fuel mixture issues, and a surrounding coaxial 
square duct, from which the coflow air flows. The inner jet’s inner radius is r\ = 0.162 cm and 
its wall thickness is wjet = 0.028 cm. The square duct’s width is 7.62 cm. Details of the burner 
construction and operation are provided in [8]. Thirty 70% CH 4 flames have been studied at three 
inlet fuel velocities (v F = 23 cm/s, 46 cm/s, and 81 cm/s) and five coflow velocities (peak coflow 
velocity t ! o, max approximately equals to 18, 39, 48, 54, and 71 cm/s) under both microgravity (pg) 
and normal gravity (1 g). Since each flame is surrounded by an air coflow and takes at most 4% 
of the cross-sectional area of the burner, the square duct is approximated as a coaxial tube with 
an identical cross-sectional area (radius ro = 4.288 cm); see Fig. 1 (left). The velocity profile 
of the fuel stream is parabolic, with the average Reynolds number in the fuel jet ranging from 
21 to 76. Due to hardware imperfections, the flow field of the coflow air is not perfectly flat and 
has a bump near the inner tube; see Fig. 1 (right). To capture this nonideality, the measured inlet 
velocity distribution has been fitted by a hyperbolic tangent function [8,9]. Since temperature 
measurements near the fuel tube exit are unavailable, the inlet temperature is set to 298 K. 


3 Computational approach 

The numerical framework is similar to those in the authors’ previous works (see, for example, [10- 
13]) with the MC-Smooth vorticity-velocity formulation [14] employed. Specifically, the gaseous 
mixture is assumed to be Newtonian, and the diffusion of each species is approximated as Fickian. 
The flow’s small Mach number implies that pressure P can be approximated as being independent 
of location in the flame, and mixture density p can be directly obtained from the ideal gas law. 
To evaluate the divergence of the net radiative flux, the model includes an optically thin radiation 
submodel [15-17] with three radiating species (H 2 0, CO, and C0 2 ). The gas-phase chemistry 
employed in this research is the GRI 3.0 mechanism [18] with all nitrogen-containing species 
(except N 2 ) removed, leaving 35 species and 217 reactions. All thermodynamic, chemical, and 
transport properties are evaluated by the CHEMKIN [19, 20] and TRANSPORT [21 , 22] subroutine 
libraries, parts of which were optimized for greater speed [23]. 
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Figure 1 : Left: Schematic of the burner used in the simulations. Right: Measured (symbols) and fitted (line) velocity 
profiles of the coflow air. 


Since each flame is axisymmetric, a two-dimensional computational domain is employed, extend- 
ing radially to r = 4.288 cm and axially to z = 12.2 cm. Except for the specific values of 
velocities, temperatures, and concentrations at the inlet, the boundary conditions are similar to 
those in [11,24] and are outlined here briefly. At the bottom boundary (inlet), v r , v z , T, and Y n 
(n = 1, ..., (Vgpecieg ) are specified to values mentioned previously in Section 2, and the definition 
of c o is applied. At the top boundary (outlet), the definition of uj is applied, and axial gradients 
of all other variables vanish. Along the axis of symmetry (/“-axis), v r and u> vanish; dv z /dr, 
dT I dr, and dY n /dr (n = 1, ..., N svec i es ) are set to zero. The boundary conditions on v z , T, and 
Y n are discretized as discussed in [11], in order to avoid convergence difficulties and/or discretiza- 
tion inaccuracies. At the outer radial boundary, v r , v z , and dY n /dr (n = 1, ..., N spec]es ) vanish, 
T = 298 K, and the definition of u is applied. 

The grid on the axisymmetric domain is a nonuniform tensor product mesh containing 129 points 
in the r direction and 202 points in the z direction. The grid points are clustered towards the burner 
surface and the centerline in order to capture the sharp gradients in these regions. The governing 
equations and boundary conditions were then discretized on this tensor product grid via standard 
nine-point finite difference stencils, which are second-order accurate in the regions of the grid that 
are equispaced and are between first- and second-order accurate elsewhere. The resulting set of 
fully coupled, highly nonlinear equations is then solved simultaneously at all grid points using a 
damped modified Newton’s method [25, 26] with a nested Bi-CGSTAB linear algebra solver [27]. 
To aid in convergence of Newton’s method, pseudo-transient continuation is performed until the 
adaptively chosen pseudo-time step exceeds a cut-off value, after which the steady-state equations 
are finally solved to a Newton tolerance of 10 -4 . All calculations were performed on workstations 
with 3.0-GHz Intel Xeon E5472 processors, and the typical memory usage for a computation on 
the 129 x 202 grid is around 2.2 GB of RAM. 
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Figure 2: Typical results of computed and measured luminous flame shapes illustrating effects of inlet fuel velocity. 
Color contours are measured luminous flame images, and white isopleths are computed luminous flame boundaries 
defined by contours of 1% of the maximum CH concentration. Flame parameters are shown at the top of each plot. 


4 Experimental approach 

A digital single lens reflex camera was used for flame shape measurement and was fully char- 
acterized as a ratio pyrometer for soot temperature measurement. The spectral responses of the 
red, green, and blue channels have been measured following the approach detailed in [28, 29]. A 
calculated lookup table that correlates color ratio and soot temperature is used to determine the 
axisymmetric soot temperature distribution. Further details can be found in [8]. 


5 Preliminary results and discussion 

The numerical model has been validated for buoyant and nonbuoyant coflow laminar diffusion 
flames in previous research [8, 9]. In this work, additional validation was conducted by comparing 
the computed and measured luminous flame shapes of a series of CH 4 -air coflow laminar diffusion 
flames in both /i g and 1 g. Figure 2 shows a comparison between computed and measured lumi- 
nous flame shapes at identical fuel dilution level to demonstrate the effects of inlet velocity and 
gravity on flame geometry. This direct comparison was enabled by the previous observation [2] 
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Figure 3: Left: Effect of inlet fuel velocity and inlet coflow velocity on modified flame length. Right: Effect of inlet 
fuel velocity and inlet coflow velocity on flame radii. In both plots, closed symbols denote measurements, and solid 
lines with open symbols denote predictions; squares denote results at 1 g, and diamonds denote results at /./, g. Plots 
(a), (b), (c), (d), and (e) correspond to data with peak coflow velocity set to 18, 39, 48, 54, and 71 cm/s. 


that excellent structural agreement exists between computed CH and measured CH* in CHzt-air 
diffusion flames for flame shape and lift-off height in both /i g and 1 g. Since the flames studied 
in our work range from lifted, non-sooting flames to attached, sooting flames, we defined the com- 
puted luminous flame boundary as the contour of 1% of the maximum CH concentration to ensure 
data consistency. Regardless of the presence of soot, there is generally good agreement between 
computed and measured flame shapes: the computed luminous flame boundaries are located very 
close to the measured ones, and they capture the tendency that the flames become longer and wider 
as buoyancy is eliminated or fuel concentration is increased. 

Figure 3 shows the effects of inlet fuel velocity, inlet coflow velocity, and gravity on the modified 
flame length (obtained by subtracting the lift-off height from the flame length) and flame radius 
of coflow laminar methane-air diffusion flames. From Fig. 3, we can see that in most flames, the 
agreement between computations and experiments is very good. In addition, we can also find that 
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Figure 4: Effects of gravitational acceleration on (a) modified flame length and (b) flame radius of 70% CH 4 flames. 
In both plots, closed symbols denote measurements, solid lines with open symbols denote predictions. Results with 
inlet fuel velocities equal to 23 cm/s, 45 cm/s, and 81 cm/s are denoted by squares, diamonds, and circles. 


for a fixed coflow velocity, increasing the fuel velocity will effectively increase the length and 
the radius of the flame. On the other hand, increasing the coflow velocity will reduce both the 
modified flame length and the flame radius. Furthermore, it can be seen from Fig. 3 that changing 
gravity also has a significant impact on laminar diffusion flames. Specifically, when gravitational 
acceleration is eliminated, the flame becomes longer and wider, and the actual difference (i.e., on 
an absolute basis) between n g flame and 1 g flame will increase when the inlet fuel velocity is 
increased or the inlet coflow velocity is reduced. 

Finally, we will take a closer look at the effects of gravity on flame structure. In Fig. 4, we 
have employed a total of five gravitational acceleration levels, 0 g, 0.02 g (typical value of the g- 
jitter on aircraft flying parabolic trajectories), 0.17 g (gravitational acceleration on Moon), 0.38 g 
(gravitational acceleration on Mars), and 1.0 g (gravitational acceleration on earth). From the 
results in Fig. 4, we see that the effects of gravity on flame structure are highly nonlinear, and 
the flame shape changes much faster near 0 g than near 1 g. In addition, we also see that this 
nonlinearity increases with increasing inlet fuel velocity, and this effect is particularly obvious in 
the flame radii results. This result is not surprising because the buoyant acceleration is countered 
by buoyant entrainment such that gravity tends to have a limited effect on the flame height. In 
contrast, buoyant entrainment tends to narrow these flames. 

Further numerical investigations are in progress. In particular, the effect of fuel stream dilution 
(e.g., 40% CH 4 and 100 % CH 4 ) will be explored at a series of coflow velocities by comparison 
with corresponding experimental measurements. After that, the sooting behavior of these flames 
will also be examined. 
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